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We report on searches for B~ Df(j> and B~ —> Dl~(f>. In the context of the Standard Model, 
these decays are expected to be highly suppressed since they proceed through annihilation of the b 
and u quarks in the B~ meson. Our results are based on 234 million T(45') ^ BB decays collected 
with the BABAR detector at SLAG. We find no evidence for these decays, and we set Bayesian 
90% confidence level upper limits on the branching fractions B{B~ —> Dftf) < 1.9 x 10“® and 
B(B~ D*~ tf) < 1.2 X 10“®. These results are consistent with Standard Model expectations. 


PACS numbers: 13.25.Hw, 12.60.Jv, ll.30.Pb 

In the Standard Model (SM), the decay B~ (p 

occurs through annihilation of the two quarks in the 
B meson into a virtual W as shown in Figure ^ No 
pure kF-annihilation B decays have ever been observed. 
The current upper limits on the branching fractions of 
B- Djcj) and ^ are 3.2 x IQ-'i (90% C.L.) 
and 4 x 10“"* (90% C.L.), respectively, set by the CLEO 
collaboration in 1993 P|. 



In the SM, B annihilation amplitudes are highly sup¬ 
pressed. Calculations of the B~ ^ Djp branching frac¬ 
tion give predictions of 3 x 10“^ using a perturbative 
QCD approach 0, 1.9 x 10“® using factorization Q, and 

7 X 10“^ using QCD-improved factorization Q- 

Since the current experimental limits are about three 
orders of magnitude higher than the SM expectations, 
searches for B ^ Dg 4> could be sensitive to new 
physics contributions. Reference 0 argues that the 
branching fraction for B~ D~(j) could be as high as 

8 X 10“® in a two-Higgs doublet model and 3 x 10“'^ in the 
minimal supersymmetric model with i?-parity violation, 
depending on the details of the new physics parameters. 

Our results are based on 234 x 10® T(4S') ^ BB de¬ 
cays, corresponding to an integrated luminosity of 212 
fb“^, collected between 1999 and 2004 with the BABAR 
detector 0 at the PEP-II B Factory 0 at the Stan¬ 
ford Linear Accelerator Center. A 12 fb“^ off-resonance 
data sample, with a center-of-mass (CM) energy 40 MeV 
below the T’(4S') resonance peak, is used to study contin¬ 


uum events, e+e” ^ qq (q = u, d, s, or c). The number of 
B mesons in our data sample is two orders of magnitudes 
larger than in the previously-published search [ij. 

We search for the decay B~ Dg </) in the fol¬ 
lowing modes: D*~ D~j, Dj —> (i)'K~, KgK~, and 
K*^K~ , with the secondary decay modes p —> K^K~, 
Kg 7r+7r“, and K*^ AT+tt”. (Charge-conjugate 
decay modes are implied throughout this article.) We 
denote the p from B ^ Dg p as the “bachelor p” in 
order to distinguish it from the p in D~ —>■ pTT~. 

Our analysis strategy is the following. Starting from 
the set of reconstructed tracks and electromagnetic 
calorimeter clusters, we select events which are kine¬ 
matically consistent with the T(4S') ^ BB, B~ 

Dg p hypothesis. Backgrounds, mostly from contin¬ 
uum events, are suppressed using a likelihood constructed 
from a number of kinematical and event shape variables. 
For each candidate satisfying all selection criteria, we cal¬ 
culate the energy-substituted mass, ruES, defined later in 
this article. The mES distribution of these events is then 
fit to a signal plus background hypothesis to extract the 
final signal yield. 

All kaon candidate tracks in the reconstructed decay 
chains must satisfy a set of loose kaon identihcation cri¬ 
teria based on the response of the internally-reflecting 
ring-imaging Cherenkov radiation detector and the ion¬ 
ization measurements in the drift chamber and the silicon 
vertex tracker. The kaon selection efficiency is a func¬ 
tion of momentum and polar angle, and is typically 95%. 
These requirements provide a rejection factor of order 
10 against pion backgrounds. No particle identification 
requirements are imposed on pion candidate tracks. 

We select p, Kg, and AT*® candidates from pairs of 
oppositely-charged tracks with invariant masses consis¬ 
tent with the parent particle decay hypothesis and con¬ 
sistent with originating from a common vertex. The in¬ 
variant mass requirements are ±10 MeV (^ 2.4r) for the 
p, ±9 MeV (- 3cr) for the AT®, and ±75 MeV (- 1.5r) 
for the A'*®, where a and F are the experimental and 
natural width, respectively, of these particles. (Here, and 
throughout the paper, we use natural units where c = 1). 
We then form Dj candidates in the three modes listed 
above by combining p. Kg, or AT*® candidates with an 
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additional track. The invariant mass of the D~ candi¬ 
date must be within 15 MeV (~ 3cr) of the known Dj 
mass. In the D~ (j)'K~ and D~ K*^K~ modes, 
all three charged tracks are required to originate from a 
common vertex. In the D~ K^K~ mode, the /fg and 
D~ vertices are required to be separated by at least 3 
mm. This last requirement is very effective in rejecting 
combinatorial background and is 94% efficient for sig¬ 
nal. We select candidates from Dj and photon 
candidates. The photon candidates are constructed from 
calorimeter clusters with lateral profiles consistent with 
photon showers and with energy above 60 MeV in the 
laboratory frame. We require that the mass difference 
AM = M{D*~) —AI{D~) be between 130 and 156 MeV. 
The AM resolution is about 5 MeV. 

At each stage in the reconstruction, the measurement 
of the momentum vector of an intermediate particle is 
improved by refitting the momenta of the decay prod¬ 
ucts with kinematical constraints. These constraints are 
based on the known mass 0 of the intermediate particle 
and on the fact that the decay products must originate 
from a common point in space. 

Finally, we select B candidates by combining Dg 
and bachelor (f) candidates. A B~ candidate is char¬ 
acterized kinematically by the energy-substituted mass 

Toes = Po- pbY/EI- and energy differ¬ 

ence AA = Ag — 5 Vi, where E and p are energy and 
momentum, the asterisk denotes the CM frame, the sub¬ 
scripts 0 and B refer to the initial T(4S') and B candi¬ 
date, respectively, and s is the square of the CM energy. 

In the CM frame, tues reduces to mss = \l— p*^. 
For signal events we expect tues ~ Mb, the known B~ 
mass, and /S.E ~ 0. The resolutions of tues and AA are 
approximately 2.6 MeV and 10 MeV, respectively. For 
events with more than one B~ candidate, we retain the 
candidate with the lowest computed from the mea¬ 
sured values, known values, and resolutions for the Dj 
mass, the bachelor (p mass, and, where applicable, AM. 

This analysis was performed blind: GEANT4 simu¬ 
lated data 0 or data samples outside the fit region were 
used for background studies and selection criteria opti- 
mization. Most of the backgrounds to the B~ ^ Dg ' cp 
signal were determined to be from continuum events. To 
reduce these backgrounds we make two additional re¬ 
quirements. First, we require |cos0t| < 0.9, where 9t 
is the angle between the thrust axis of the B~ candi¬ 
date and the rest of the tracks and neutral clusters in 
the event, calculated in the CM frame. The distribution 
of I cos 6t I is essentially uniform for signal events and 
strongly peaked near unity for continuum events. Sec¬ 
ond, for each event we define a relative likelihood for 
signal and background based on a number of kinematical 
quantities. The relative likelihood is defined as the ratio 
of the likelihoods for signal and background. The sig¬ 


nal (background) likelihood is defined as the product of 
the probability density functions, PDFs, for the various 
kinematical quantities in signal (background) events. 

The kinematical quantities used in the likelihood are 
reconstructed masses, helicity angles, and a Fisher dis¬ 
criminant designed to distinguish between continuum 
and BB events; we will discuss each of these in turn 
in the following paragraphs. All PDFs are chosen based 
on studies of Monte Carlo and off-resonance data. 

The masses used in the likelihoods are those of the 
Dj, AM for D*- Dg-f, the K*° in Dj K*°K-, 
and the (p in Dj (pTv~. The signal PDFs for the mass 
variables are the sum of two Gaussian distributions for 
Dj and AM, a Breit-Wigner distribution for the K*^, 
and a Voigtian distribution for the p. We parameter¬ 
ize the background PDFs as uniform distributions. Note 
that the mass of the bachelor p and the mass of the 
in Dj K°K- are not used in the definition of the 
likelihoods. This is because studies of background event 
samples suggest that background events contain mostly 
real bachelor p and real mesons. 

The helicity angles used in the likelihood are those in 
K*o iP+TT" and in ^ > K^K~ for both p candidates. 

The signal PDFs for these quantities are set by angu¬ 
lar momentum conservation to be proportional to cos^ 6, 
where 9 is the helicity angle for the process. The one ex¬ 
ception is the helicity angle distribution of the bachelor p 
in B~ —> D*~p, where the polarization of the two vector 
mesons in the final state is not known. For this reason, 
the helicity angle of the bachelor p is not used in the 
definition of the likelihood for the B~ D*~p mode. 
The background PDFs for these variables are uniform in 
cos 6*. In addition, in the likelihood we also use the polar 
angle of the B~ candidate in the CM frame {9b)- The 
signal is expected to follow a sin^^s distribution, while 
the background is independent of cos 9b- 

The final component of the likelihood is a Fisher dis¬ 
criminant constructed from the quantities Lq = J^iPi 
and L 2 = Y^^p* cos^ a*. Here, p* is the magnitude of 
the momentum and a* is the angle with respect to the 
thrust axis of the B~ candidate of tracks and clusters not 
used to reconstruct the B~, all in the CM frame. The sig¬ 
nal and background PDFs for this variable are modeled 
as bifurcated Gaussians with different means and stan¬ 
dard deviations. Note that the likelihood is constructed 
to have negligible correlations between its components. 
Since the Fisher discriminant is highly correlated with 
the I cos 9t I variable defined above, the | cos 9t \ variable 
is treated separately and not included in the likelihood. 

The likelihood requirements have been optimized to 
maximize the expected sensitivity on a mode-by-mode 
basis. (The expected sensitivity is defined as the branch¬ 
ing ratio upper limit that we should obtain, on average, 
when performing this experiment many times in the ab¬ 
sence of a signal.) Depending on the mode, the combined 
efficiency of the likelihood and the | cos 9t \ requirements 
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TABLE I: Efficiencies (e), branching fractions {B), and prod¬ 
ucts of efficiency and branching fractions for the modes used 
in the B~ —> search. The uncertainties on the e and B 

are discussed in the text. Here B is the product of branching 
fractions for the secondary and tertiary decays in the specified 
decay mode. 


B Mode Dg 

Mode 

e 

H(10-'’) 

e X H(10-'’) 

B ^ Dg(f> 

D7 ^ 

(pTV 

0.192 

11.6 

2.22 

D7 ^ 

K-Kl 

0.177 

8.2 

1.45 

D7 ^ 

K*°K- 

0.140 

14.5 

2.03 

B- ^ orp 

D7 ^ 

(pTV 

0.109 

10.9 

1.19 

D7 ^ 

K-K° 

0.100 

7.7 

0.77 

D7 ^ 

K*°K- 

0.083 

13.6 

1.14 


varies between 71% and 83%, while providing a rejec¬ 
tion factor of between 4 and 7 against continuum back¬ 
grounds. 

After applying the requirements on relative likelihood 
and I cos 9t\, we also demand that AE fall inside the sig¬ 
nal region: within 30 MeV 3cr) of its expected mean 
value for signal events. This mean value is determined 
from simulation, and varies between —3 and OMeV, de¬ 
pending on the mode. 

The efficiencies of our selection requirements, shown 
in Table Q] are determined from simulations. For the 
B~ D*~(j) mode, we take the average of the efficien¬ 
cies calculated assuming fully longitudinal or transverse 
polarization for the two vector meson final state. These 
efficiencies are found to be the same to within 1%. The 
quantities B in Table Q] are the product of the known 
branching fractions for the secondary and tertiary decay 
modes. These are taken from the compilation of the Par¬ 
ticle Data Group 0 , with the exception of the branching 
fraction for D~ for which we use the latest, most 

precise measurement B{D~ (/'tt”) = (4.8 ± 0.6)% [3. 
Since the branching fractions for the other two Dj modes 
are measured with respect to the Dj (f>Tr~ mode, we 
have rescaled their tabulated values from the Particle 
Data Group accordingly. 

The systematic uncertainties on the products of effi¬ 
ciency and branching ratio for the secondary decays in 
the decay chain of interest are summarized in Table m 
The largest systematic uncertainty is associated with the 
uncertainty on the D~ —>■ 4n:~ branching ratio, which is 
only known to 13% |9|, and which is used to normalize 
all other Dj branching ratios. 

The experimental systematic uncertainties all relate to 
the determination of the efficiency, e. The dominant 
source of error is the uncertainty in the efficiency of 
the kaon identification requirements. The efficiency of 
these requirements is calibrated using a sample of kine¬ 
matically identihed ^ iD°7r+, K~Tr'^ de- 


TABLE II: Systematic uncertainties on Ci ■ Bi, where the 
index i runs over the three D~ modes used in this analysis, 
ei are the experimental efficiencies, and Bi are the branching 
fractions for the i**' mode. 


Source 

B —> Dg (j) B 

- ^ Drp 

D 7 branching fraction 

14% 

14% 

D*~ branching fraction 

- 

2.5% 

Other branching fractions 

1.5% 

1.5% 

Charged kaon ID 

13.2% 

13.3% 

Tracking and efficiency 

3.7% 

3.7% 

Photon efficiency 

- 

1.8% 

Final state polarization 

- 

1% 

Selection requirements 

5% 

5% 

Simulation statistics 

0.6% 

0.6% 

Total 

20% 

21% 


cays, where track quality selection differences between 
this sample and our analysis sample have been taken 
into account. The efficiency of the kaon identihcation re¬ 
quirements is assigned a 3.6% systematic error. This re¬ 
sults in a systematic uncertainty of 14% for the efficiency 
of the modes with four charged kaons {D~ K*^K~, 

D~ and 11% for the mode with three charged 

kaons {Dj K^K~). A second class of uncertainties 
is associated with the detection efficiency for tracks and 
clusters. From studies of a variety of control samples, 
the tracking efficiency is understood at the level of 1.4% 
(0.6%) for transverse momenta below (above) 200 MeV. 
There is also a 1.9% uncertainty associated wih the re¬ 
construction of iFg —*■ TT+TT" which can occur a few cen¬ 
timeters away from the interaction point. Given the mul¬ 
tiplicity and momentum spectrum of tracks in the decay 
modes of interest, the uncertainty on the efficiency of re¬ 
constructing tracks in the i?-decay chain is estimated to 
be 3.7%. In the B~ search, there is an addi¬ 

tional uncertainty of 1.8% due to the uncertainty on the 
efficiency to reconstruct the photon in D*~ —> Dgy, and 
a 1% uncertainty from the unknown polarization in the 
final state. Finally, to ascertain the systematic uncer¬ 
tainty due to the efficiency of the other event selection 
requirements, we examine the variation of the efficiency 
under differing conditions: shifting the AE by 3 MeV 
(0.3%); shifting the mean of the Dj and (j) masses and 
AM by 1 MeV (0.2%, 0.1%, 0.2%, respectively); in¬ 
creasing the width of the Dj and (p masses and AM by 
1 MeV (1.5%, 0.4%, 1.5%, respectively); using a Fisher 
distribution obtained from the data sample of a similar 
analysis, B Dtt with D Ktt (3%). Thus we as¬ 
sign a 5% systematic on the combined efficiency of these 
selection criteria. 

We determine the yield of signal events from an un¬ 
binned extended maximum-likelihood fit to the toes 
distribution of B~ candidates satisfying all of the re¬ 
quirements listed above. We ht simultaneously in two 














7 


> 


Qi 

§ 


VI 


> 





mEs(GeV) 



mEs(GeV) 


FIG. 2: Distribution of rriES for (a) B” —> D~(l) and (b) 
B~ —> Dl~(j) candidates in the \AE\ signal region. The su¬ 
perimposed curves are the result of the fits described in the 
text. The dashed curve is the background contribution and 
the solid curve is the sum of the signal and background com¬ 
ponents. 


|Ai?| regions: in the signal region the distribution is 
parametrized as a Gaussian and the combinatorial back¬ 
ground as a threshold function na ; in a sideband of AE 
(|Ai?| < 200 MeV, excluding the signal region) we fit 
solely for the threshold function parameter. In our fit, the 
amplitude of the Gaussian is allowed to fluctuate to neg¬ 
ative values, but, for reasons of numerical stability, the 
sum of the Gaussian and the threshold function is con¬ 
strained to be positive over the full tties fit range. The 
mean and the standard deviations of the Gaussian are 
constrained to the values determined from Monte Garlo 
simulation. The fitting procedure was extensively tested 
with sets of simulated data, and was found to provide an 
unbiased estimate of the signal yield. 

Figure El shows the mss distribution of the selected 
candidates. We see no evidence for B Dg (j)- Tfi® 
fitted event yields are N = —I.G^qq and N = for 

the B~ D~<j) and B~ modes, respectively, 

where the quoted uncertainties correspond to changes of 
0.5 in the log-likelihood for the fit. The likelihood curves 
are shown in Figure |21 The requirement that the sum 
of the Gaussian and the threshold function be always 
positive results in an effective constraint N > —1.6 in the 
B~ D~(j) mode. This is the source of the sharp edge 
at TV = —1.6 in the likelihood distribution of FigureEJa). 

We use a Bayesian approach with a flat prior to set 90% 
confidence level upper limits on the branching fractions 
for the B~ D~(j> and B~ —s- D*~(p modes. In a given 
mode, the upper limit on the number of observed events 



FIG. 3: Likelihood from the fit in arbitrary units as a function 
of the number of signal events, (a) B~ Dj(j>\ (b) B~ 


{Nul) is defined as 


rNuL 


C{N) dN = 0.9 




C{N) dN. (1) 


C{N) is the likelihood, determined from the mss fit de¬ 
tailed above, as a function of the number of signal events, 
N. The upper limit B on the branching fraction is 


B < 


Nul 

^BB Si ^ 


( 2 ) 


^BB ~ (233.9 ± 2.5) X 10® is the number of BB events, 
the index i runs over the three D~ decay modes, is the 
efficiency in the mode, and Bi is the product of all 
secondary and tertiary branching fractions (see Table QJ. 

We account for systematic uncertainties by numeri¬ 
cally convolving C{N) with a Gaussian distribution with 
its width determined by the total systematic uncertain¬ 
ties (Table HH) in the two modes, including the 1.1% un¬ 
certainty in added in quadrature. We find limits 

B{B- D-(j)) < 1.9 X 10-® and B{B- < 

1.2 X 10“® at the 90% confidence level. These limits are 
calculated using B{Dj = (4.8 ± 0.6)% from 

Reference |i|. If we were to use the value B{D~ 
(j)TT~) = (3.6 ± 0.9)% from the Particle Data Group |^, 
we would find B{B~ D~(j)) < 2.7 x 10“® and B{B~ 

< 1.7 X 10“®. For completeness, we also com¬ 
pute B{B~ Djcj)) X B{Dj — > ^7r“) < 8.6 x 10“® and 
B{B~ —> D*~(j}) X B{Dj < 5.4 x 10“^, also at 

the 90% confidence level. 

In summary, we have searched for B —>■ Dg ' (f>, and 
we have found no evidence for these decays. Our limits 
are about two orders of magnitude lower than the previ¬ 
ous results, but are still one order of magnitude higher 
than the Standard Model expectation. 
































































Using the calculation of Reference 0 , we can use our 
results to set bounds on new physics contributions. These 
bounds are obtained in the framework of factorization, 
neglecting systematic uncertainties associated with the 
calculation of hadronic effects. In the type II Two- 
Higgs-Doublet model we extract a tree-level 90% C.L. 
limit tan/3/MH < 0.37/GeV, where tan/3 is the ratio of 
the vacuum expectation values for the two Higgs dou¬ 
blets and Mh is the mass of the charged Higgs boson. 
This limit is not quite as stringent as the limit that 
can be obtained from the B~ t~Vt decay mode, 
tan/3/MH < 0.29/GeV M 

In the context of supersymmetric models with i?-parity 
violation (RPV) the new physics contribution to 

B ^ Dscj) depends on the quantity , 

where is the coupling between the generation 

doublet lepton superfield, the generation doublet 
quark superfield, and the generation singlet down- 
type quark superfield; Mi is the mass of the generation 
charged super-lepton. Conservatively assuming maximal 
destructive interference between the SM and RPV am¬ 
plitudes, we find |^| < 4 x 10“'^/(100 GeV)^. 
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